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1. Introduction 

Few-body hadron-hadron interaction has been and is still a fundemantal constituent part of 
hadron and nuclear physics. Among them, the antinucleon-nucleon (NN) has been achieved fruit¬ 
ful progress, especially in a meson-exchange model, see e.g., a review in Ref. ||I]]. And after 1990s, 
chiral effective field theory (EFT) has become a powerful tool to analyze nucleon-nucleon interac¬ 
tion, for a review, see e.g.. Ref. [Q]. However, only little work of NN in chiral EFT has been done, 
besides the one of partial-wave analysis (PWA) |^. The recently resurgent interest of pp physics 
is triggered by the threshold enhancement in pp invariant mass spectrum observed in experiments, 
e.g, for the decays J/y —> Ypp 0, ||], pp |Q, ^]. We will elaborate below how we use the 

pp final-state interaction (FSI) to interpret these phenomenons. 

2. NN interaction in chiral EFT 

NN potential consists of two parts: elastic scattering part and the annihilation, which is the 
same feature in the framework of both the conventional meson-exchange model and the chiral 
EFT. The difference is on the technical treatment. In chiral EFT, the elastic part is governed by 
the pion exchanges (pion as the only degree of freedom in chiral EFT), which is tied closely to the 
knowledge of NN interaction except for the sign difference due to the G—parity transformation. 
For the power counting rule and details in NN scattering, one refers to Refs. 0. Here we only 
remind that 


yNN ^ _yNN yNN ^ yNN 


(2.1) 


because of the G—parity transformation rule from nNN vertex to nNN. The new feature of NN 
compared to NN is the existence of annihilation effect that will be parameterised in contact term (in 
charge of short-range interaction). In liilich model 0, the annihilation is treated as a energy-, spin- 
, and isospin-independent Gaussian form. Here we still follow the spirit of chiral power counting, 
taking partial wave as example. Starting form the most elaborated couple-channel model, we 
have 


Knn — ^ ^NN^X^x{z)Vx^nN^ 

X=2n,'in.- 

where X denotes, in principle, any possible intermediate states including 2?!, 3?!, etc., and G is 
the free Green’s function. It is argued that the annihilation does not introduce a new scale into the 
problem |^], i.e., it can be likewise treated in the chiral expansion. Then up to next-to-next-to- 
leading order (NNEO), one can write the NN i (i as mesons) annihilation potential in 'Sq partial 
wave as 

= + (2-3) 

where p{p') is the module of three-momentum in center-of-mass system (CMS) of intital (final) 
NN states. Picking out the imaginary part, one would get 

ImV^nCSo) = - (cf,^+a,y) . (2.4) 
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Equation (2.4) fulfils the unitary condition by definition (an alternative method is based on the 
dispersion theory, see e.g., the nn sector in Ref. [[I^]). Expanding the real part from the principal 
integral, we will get the similar structure as the NN case. 


ReKnn('5o) = +C,5„ {p^ + p'^) . (2.5) 

Note there are four EECs (Cf^^, Ci^^, Ci^^) in total, where the annihilation in indicated in 

the superscript by “a”. Now the NN potential is setting up but containing the low-energy con¬ 
stants (EECs). The recent energy-dependent NN partial-wave analysis (PWA) is done in Ref. |Q], 
which provides a rather nice description of all the pp scattering data below laboratory momen¬ 
tum of 925 MeV. These EECs will be fitted to the partial-wave amplitude there. The results for 
the inelasticity and phase shfits of are shown in Eig. H up to Tiab = 200 MeV for NEO and 
Tiab = 250MeV for NNEO. Tiab is the kinetic laboratory energy, Tiab = 2k^/m with k denoting 
the module of the on-shell momentum in CMS. We have used the notation where 

L, S, J denote the orbital angular momentum, total spin and their quantum addition, i.e, total angu¬ 
lar momentum, respectively. The phase shift (complex value due to annihilation) is defined from 
5—matrix as S = with 5 = 5r + i5j, and then 5/ = — log(T])/2. The band is formed by 

varying the cutoff combination applied into the Eippmann-Schwinger equation and spectral func¬ 
tion in the two-pion exchange potential |[^, Erom Eig. [Tjone sees that our results reproduce 
the PWA rather well with very small uncertainty (cutoff dependence) at the whole region consid¬ 
ered. The corresponding results for coupled partial wave — ^D\ will be used in the following 
subsection but are not shown here due to the limited room. Eor the reader who is interested in this 
part, one refers to the publication in [^]. Besides the phase shifts and inelasticities, the scattering 
length, and the level shifts and widths of the antiprotonic hydrogen are also calculated. They are 
all in good agreement with experimental numbers within error bars. NN Bound states are predicted 
in and partial waves 0. 


3. pp related hadron physics 

3.1 pp-threshold enhancement in 7/i/r —> ypp 

After the discovery of the strong threshold enhancement observed in //i/r —)• ypp by BES 
collaboration 0, many explanations have been proposed. Due to its proximity to pp threshold, it 
is speculated to be a pp bound state, or at least, has much to do with pp interaction. To take into 
account the pp ESI, we write the total amplitude symbolically as 


A —Aq +AqGTpp, 


(3.1) 


where Aq is the elementary production amplitude without considering pp ESI, and in the second 
term the off-shell form for Aq is needed since it appears in the integral; G is the free pp Green’s 
function; the P-matrix elements, Tpp, can be calculated from Sec. Writing Eq. ( pA| ) in a partial 
wave more explicitly, one gets (here for S—wave) 


Al=aI 



dqq^ 1 

( 271)3 lEk-2Eq + iQ+ 


TL{q,k-,Ek) 


(3.2) 
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Figure 1: Complex phase shifts for isospin-0 '5o (^*5o) partial wave in unit of degree as a function of 
kinetic laboratory energy Tiab- 5^ coincides with the widely used conventional real phase shift while 5j = 
— log(T7)/2. Circle points represent PWA reported from Ref. [^. Green (Red) band indicates the cutoff 
dependence of NLO (NNLO) potential. 


At low-energy region, assuming Aq has only a weak energy dependence, one may reasonably ap¬ 
proximate as a constant. In this way, in the model there will be only one parameter, i.e., the overall 
normalization constant. For the observable of event distribution, we are doing a parameter-free 
calculation, as a matter of fact, in viewpoint of only the energy dependence. We stress that the 
energy dependence of this whole system solely comes from pp interaction since Aq is treated as 
a constant. This point indeed verifies the momentous role of pp FSI. In the channel pp, 

see Sec. |3.2[ , the overall constant is used to match the magnitude of the cross section. In the pro¬ 
cess y/i/r —)■ ypp, the lowest allowed quantum number for pp is 'Sq, and both isospin-0 and 1 are 
allowed. The result based on the original ^Sq potential constrained by PWA of Ref. does not re¬ 
produce such a strong enhanced peak near pp threshold. Instead, we perform a combined analysis 
of the pp scattering data as well as the prominent peak shown by BES data [^. T-matrix is taken 
as r = (r° + r^)/2, where the superscript denotes isospin. The LECs for isospin-0 is kept 
as what comes from fitting to the PWA of Ref. [Q], supported by the milder energy dependence of 
y/t/r —)• (opp [Vl\, while the four LECs in isospin-1 is refitted. So in total, there are 5 free pa¬ 
rameters (4 LECs -I- 1 overall constant). The results are presented in Eig. ^and Eig. One can see 
that the peak around pp threshold is nicely reproduced, and simultaneously, the 'Sq partial-wave 
cross section of the PWA as well the original one constructed by us but fitted to PWA are very 
well reproduced. Then the total cross section from our such potential (only part changes and 
others do not alter compared to Ref. [^) is thus expected to be in a good agreement with the one 
calculated from PWA. The protonium level shift and widths are also examined, and they are also 
within experimental error bars. It turns that in order to describe a such prominent peak, we need a 
pp bound state in isospin-1 However, the data above threshold is believed to be not capable of 

pinning down the binding energy and width of this bound state. More information on the invariant 
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Figure 2: pp spectrum for the decay J/y/ ^ ypP- The band represents our result. The dashed curve denotes 
the phase space behavior. Data are taken from Refs. ^ [^]. The measurement of Ref. [|[] is adopted for 
the scale. The data for the BES measurement from 2003 have been shifted 1 MeV to the right to discriminate 
from the new measurement. 


mass spectrum below pp threshold is needed. In Ref. [[I4|], the systematic description of the pp 
mass spectrum in other J/y and y' decays are achieved. In our calculation, the mass difference 
of proton and neutron, as well as the Coulomb interaction is not considered, and further work is 
ongoing. 


3.2 Low-energy e^e —)• pp observables 


Examining electromagnetic form factors (EMFF) of the proton {Ge and Gm) is an efficient 
way to probe the nucleon structure. The reaction PP, and its inverse one pp ^ 

(these two are related to each other by time reversal operation) are used to measure EMFF. The 
experiment shows a strong energy dependence in proton EMFF close to pp threshold. Recent 
measurements were done in Refs. ij- As shown above, at such energy region, the pp FSI plays 
an important role. And here in the > pp decay, it is no exception. Taking into the fact 

that one photon exchange should dominate in the decay pp, and thus the only allowed 

partial wave is the coupled — ^Di. This provides an opportunity to make a (somewhat) clean 
prediction, compared to other decay channels where many partial waves are possible and maybe 
have a comparable significance. Including the coupled partial wave — ^D\, one could extend 
Eq. (0) to a 2 X 2 matrix form 


Aq 


(Af Af \ 



(3.3) 


The matrix Aq, again as before, is the bare production amplitude without pp FSI, and is connected 
to the bare EMFF G^ and G^. At near-threshold region, they can be approximated as constant. And 
imposing the condition Ge = Gm we have only one overall normalization constant. Concerning the 
NN potential used in this work, both the chiral potential constructed in Ref. [Q] and Jiilich model 
A(OBE) are considered. The results are shown in Fig. ^ for the cross section and effective 
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Figure 3: The partial-wave cross sections as a function of the excess energy. The squares represent the 
results for the published NNLO potential 0 with the cutoff combination {450 MeV, 500 Me V}. The circles 
indicate the cross sections for the partial-wave amplitudes of Ref. [|^. The bands show the results based on 
the refitted isospin-1 *5o amplitudes. 


form factor, where for the cross section we fitted to 60 MeV, and the effective form factor are 
calculated from the fitted overall constant. The ratio and the phase difference between the proton 
form factors Ge and Gm are also presented up to the same energy region as the cross section, see 
Ref. [|I7|]. As can be seen, it reproduces the data rather well. We also calculate the differential cross 
section at a lower excess energy of 36.5 MeV, and the data is nicely reproduced, see Fig. ||. These 
altogether provide a good description of low-energy data on -^PP, and thus strongly support 
our speculations of large pp FSI. In the reactions ^multipions, pp is also shown to play an 
important role in the region [1750, 1950] MeV [18|. 


4. Summary 

In summary, we have constructed and established a NN potential in chiral EFT. The resulting 
phase shifts and inelasticities agree with the partial-wave analysis reported in Ref. |j^. Scattering 
lengths and the level shift and widths of antiprotonic hydrogen are calculated, and they are in a good 
agreement with the experimental information within error bars. With such a NN interaction at hand, 
we explored the pp FSI in several reactions. For J/\p ^ ypp perform a combined analysis 
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Figure 4; Cross section and effective form factor of the reaction pp as a function of the excess 

energy. The data are from the DM1 [ [f^ (triangles), FENICE [ jl^ (squares), and BaBar [||] (empty circles), 
(filled circles) collaborations. The red/dark band shows results based on the NN amplitude of the chiral 
EFT interaction up to NNLO while the green/light band are those for NLO. The solid line is the result 
for the NN amplitude predicted by the Jiilich model A(OBE) [^. The BaBar 2006 data are shifted to higher 
excess energy by 1 MeV. 


of experimental events distribution and the pp partial-wave cross section (^Sq case). In order to 
describe the prominent peak shown in J/y/ ^ ypp, we need a bound state in But the binding 
energy can not be well determined. The large pp FSI is also verified in ^ pp. The cross 
section and effective form factor up to the excess energy of 100 MeV are nicely reproduced. In all 
these processes, in fact, the whole energy dependences with the pp invariant mass spectrum come 
solely from the pp final-state interaction, since the bare production amplitude is approximated as 
constant without energy dependence. 
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Figure 5: Differential cross section for e+e —5> pp at the excess energies of 36.5 MeV. The data are an 

average over [0, 73] MeV and are taken from Refs. ^]. Same notations as in Fig. 4. 
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